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Abstract 

In String theory realizations of inflation, the end point of inflation is of- 
ten brane-anti brane annihilation. We consider the processes of reheating 
of the Standard Model universe after brane inflation. We identify the 
channels of inflaton energy decay, cascading from tachyon annihilation 
through massive closed string loops, KK modes, and brane displacement 
moduli to the lighter standard model particles. Cosmological data con- 
strains scenarios by putting stringent limits on the fraction of reheating 
energy deposited in gravitons and nonstandard sector massive relics. We 
estimate the energy deposited into various light degrees of freedom in the 
open and closed string sectors, the timing of reheating, and the reheat- 
ing temperature. Production of gravitons is significantly suppressed in 
warped inflation. However, we predict a residual gravitational radiation 
background at the level ^Iqw ~ 10~^ of the present cosmological energy 
density. We also extend our analysis to multiple throat scenarios. A viable 
reheating would be possible in a single throat or in a certain subclass of 
multiple throat scenarios of the KKLMMT type inflation model, but over- 
production of massive KK modes poses a serious problem. The problem 
is quite severe if some inner manifold comes with approximate isometrics 
(angular KK modes) or if there exists a throat of modest length other 
than the standard model throat, possibly associated with some hidden 
sector (low- lying KK modes). 
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1 Generalities: Reheating the Universe 



The transfer of inflaton energy into radiation energy in the process of (p)reheating 
after inflation is a vital part of any model of early universe inflation. According 
to the inflationary scenario, the universe at early times expands quasi-exponentially 
in a vacuum-like state without entropy or particles. During this stage of inflation, 
all energy is contained in a classical slowly moving inflaton field (f). Eventually the 
inflaton field decays and transfers all of its energy to relativistic particles, to start 
the thermal history of the hot Friedmann universe. 

The Quantum Field Theory of (p)reheating, i.e. the theory of particle creation 
from the inflaton field in an expanding universe, describes a process where the quan- 
tum effects of particle creation are not small, but instead produce a spectacular 
process where all the particles of the universe are created from the classical inflaton. 
The theory of particle creation and their subsequent thermalization after inflation 
has a record of theoretical developments within QFT. The four-dimensional QFT 
Lagrangian £(0, x, Aj, /ij^, ...) contains the inflaton part with the potential V{(j)) 
and other fields which give sub dominant contributions to gravity. In chaotic infla- 
tionary models, soon after the end of inflation the almost homogeneous inflaton field 
(j){t) coherently oscillates with a very large amplitude of the order of the Planck mass 
around the minimum of its potential. Due to the interactions of other fields with the 
inflaton in £, the inflaton field decays and transfers all of its energy to relativistic 
particles. If the creation of particles is sufficiently slow (for instance, if the inflaton 
is coupled only gravitationally to the matter fields), the decay products simultane- 
ously interact with each other and come to a state of thermal equilibrium at the 
reheating temperature T^. This gradual reheating can be treated with the perturba- 
tive theory of particle creation and thermalization p. However, for wide range of 
couplings the particle production from a coherently oscillating inflaton occurs in the 
non-perturbative regime of parametric excitation [2]. This picture, with variation in 
its details, is extended to other inflationary models. For instance, in hybrid inflation 
(including D-term inflation) inflaton decay proceeds via a tachyonic instability of the 
inhomogeneous modes which accompany symmetry breaking [S]. One consistent fea- 
ture of preheating - non-perturbative copious particle production immediately after 
inflation - is that the process occurs far away from thermal equilibrium. 



2 



Since hybrid inflation is the closest field theory prototype of string theory brane 
infiation, it will be especially instructive to refer to the theory of tachyonic preheating. 
Hybrid inflation involves multiple scalar fields $ in the infiaton sector. It can be real- 
ized in brane infiation [Hj, where the inter-brane distance is the slow rolling infiaton, 
while the subsequent dynamics of the branes are described by tachyonic instability in 
the Higgs direction. Tachyonic instability is very efficient so that the backreaction of 
rapidly generated fiuctuations does not allow homogeneous background oscillations to 
occur because all the energy of the oscillating field is transferred to the energy of long- 
wavelength scalar field fiuctuations within a single oscillation. However, this does not 
preclude the subsequent decay of the Higgs and infiaton inhomogeneities into other 
particles, and thus very fast reheating. Particles are generated in out-of-equilibrium 
states with very high occupation numbers, well outside of the perturbative regime. 

Recent developments in string theory infiation, while in the early stages, have to 
address the issue of the end point of infiation, specifically (p)reheating immediately 
after infiation. As in the QFT theory of reheating. String Theory reheating should 
be compatible with the thermal history of the universe. We can use this criterion to 
constrain the models. Yet, we are especially interested in the specific String Theory 
effects of reheating which may be, in principle, observationally testable. /^From a 
theoretical point of view, reheating after string theory inflation deals with the theory 
of particle creation in String Theory, which is an exciting topic by itself. 

In this paper we study the transfer of energy into radiation from string theory 
inflation based on brane-anti brane annihilation. Brane annihilation is a typical 
end-point of brane inflation. We assume a "warped" realization of brane inflation, 
constructed at the top of the ground of the KKLT stabilized vacuum |7] . The models 
of string theory inflation with warping branch into different possibilities. We mostly 
study reheating in the brane-antibrane warped inflation of ^3], but the methods shall 
be relevant for the models like [SI IH] • Reheating in other versions of string inflation, 
like D?) — D7 inflation of potentially can be described by QFT reheating 
while the racetrack inflation of ^2] relies on fleld theory entirely. 

The warped geometry of string theory inflation has some common features with 
the warped Randall- Sundrum flve dimensional braneworld. Based on this analogy, 
reheating in warped string theory inflation was modelled on the braneworld formalism 
in the recent paper [TH] . 
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Successful reheating means the complete conversion of inflaton energy into ther- 
mal radiation without any dangerous relics, in order to provide a thermal history of 
the universe compatible with Big Bang Nucleosynthesis (BBN), baryo/leptogenesis, 
and other observations. Dangerous relics can be massless or massive, and they are 
each a danger in their own way. Too many massless relics like radiation of gravita- 
tional waves is excluded by BBN, while too many massive relics overdose the universe. 
Therefore we have to monitor undesirable relics in string theory inflationary scenario. 
There are significant difference obetween string theory reheating and QFT reheating 
in this respect. Indeed, in string theory we expect excitation of all modes interact- 
ing gravitationally, gravitational waves, moduli fields, and KK modes, and we need 
their complete decay or extra tuning to go through the needle eye of observational 
constraints on potential non-SM particles. 

We study reheating after brane-anti-brane annihilation in warped inflation, more 
specifically, in the well-known ten dimensional model of Klebanov-Strassler (KS) 
throat geometry. We begin with a single throat case to identify systematically the 
channels of energy transfer from inflaton to radiation. There are several processes of 
energy cascading, from D — D pair annihilation and closed string loop decay to exci- 
tation of KK modes, and from them to the excitation of open string modes of the SM 
branes. Then we extend our investigation to the multiple throat case, associated with 
different energy scales. We shall estimate the decay rate of each channel, its energy 
scales, the reheating temperature and the abundance of dangerous relics. There- 
fore we have to combine the string theory picture of all relevant excitations, their 
QFT description when it is adequate, and elements of the theory of early universe 
thermodynamics. This is a challenge for string theorists and cosmologists. 

The plan of the paper is the following. 

In Section |21 we lay-out the background model, specifically KS geometry with 
extra ingredients needed for a KKLT stabilized string theory solution. We discuss 
the effect of warping on the local string modes, which will be major players in the 
reheating process. In Section IHl we study qualitatively the cascade of energy from 
D — D pairs to radiation. Section E] wraps up the thermalization process in the single 
throat model. In Section El we consider the multiple throat case. 
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2 Warped Compactification with Hierarchy 



In finding a realistic model of universe in string theory, one of more severe constraints 
come from the so-called moduli problems in cosmology. If there are light scalar fields 
around, especially those associated with moduli fields of compactification, inflation 
process could easily read to accumulation of energy in these light modes which then 
interferes with either exit from inflation or with low energy physics at later stage of 
cosmological evolution. 

A simple way out of these moduli problems, which are being investigated in string 
theory, involves turning on anti-symmetric tensor fields along compactified internal 
dimensions. In a generic situation with all possible fluxes turned on and all possible 
nonperturbative corrections included, it is beheved that the only massless degrees of 
freedom surviving the flux compactification would be that of 4-dimensional gravity 
or its supersymmetry completion in case of supersymmetric vacua. 

A common feature of these flux compactification is the warping. In this paper, we 
will be employing well-understood case of IIB compactification, where the warping 
can be summarized as the ten dimensional metric of the form, 

G = H-^l\y)g^,{x)dx^dx' + H^'\y)gijdy' dy' (2.1) 

where ^ is a Calabi-Yau metric on 6-dimensional compact manifold, and the warp 
factor H depends only on the internal Calabi-Yau direction. Note that, as far as the 
internal manifold goes, this way of writing the metric is a mere convention. Unless the 
physical process concerned depends crucially on underlying supersymmetric structure, 
we may as well rewrite the metric as 

G = H-'/\y)g^,dxV + hjjdy'dy' (2.2) 

since h = H^^^Q is the relevant physical metric of the compact direction. This form 
of the metric makes it very clear that the most important consequence of the warp 
factor, i.e., a generation of hierarchy via exponential red-shift is really coming from 
the in front of the first piece. Precise form of this warp factor is found by 

solving a second order equation for H with various source term, such as contribution 
from the NS-NS and RR fluxes, D-branes, Orientifold planes. 
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2.1 Klebanov-Strassler Throats 

An example of the warped geometry is well-understood KS solution [11] [IZI [El • KS 

geometry is one of the building block of the KKLT stabilized solution, which also 
includes warped instanton branes and Z)3 at the tip of the throat. Using this as a 
background (asymptotic of late time cosmological evolution) inflation can be realized 
by inclusion of additional elements: mobile D?> brane in the throat attracted to 
another Z)3 brane around the bottom. This — Z)3 interacting pair provides in 
inflation in 4d description, end point of inflation is their annihilation. We also need 
other D3 branes around the bottom for SM phenomenology. Geometry of the model 
is shown in Figure H 

The Klebanov-Strassler throat starts from a conifold part of the Calabi-Yau metric 
Q. Local form of the metric is conical, 

Qijdy'dy^ = dr^ + rHs^^^ (2.3) 

where the five- dimensional manifold T^'^ is a fibred over a product manifold S'^ x 
S"^. The metric is clearly singular at origin r = because the size of T^'^ vanishes 
there. 

Actual geometry involves a deformation that blows up one S*^, as a fibration of 
over one S*^, to keep it finite size at origin while allowing the other S"^ collapse to 
zero size. Thus, at the bottom r = 0, the geometry is roughly that of 

X (2.4) 

Corresponding to 5''^ is a 3-cycle on which RR 3-form flux F(3) is supported, while R^ 
is part of its dual 3-cycle on which NS-NS 3-form flux H^^,) is supported. Denoting 
these two 3-cycles as A and i?, we have the quantization condition 

with integers M and K. 

Note that the B cycle should extend outside of this local conifold geometry and the 
quantization condition on if(3) refers to this entire 3-cycle. However, assuming that 
most of if(3) flux resides within this conical region, the warp factor has been solved 
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Figure 1: Radial geometry of a Klebanov Strassler throat. For most part of this paper, 
we consider KKLMMT-hke inflation scenario where unstable D-brane system of D3 
branes and anti-D3 branes near the bottom of the throat drives inflation, possibly 
with some leftover D3's. 

explicitly. Away from the conifold point r = 0, we have the following approximate 
solution P^ . 

H=^ {rX + Rt \og{T/R^Y) (2.6) 

with 

= ^a'^g^MK, R' = A^«'2 (2.7) 

The constant part of H is to be determined by gluing this local geometry to the rest 
of the compact Calabi-Yau metric near r = But for simplicity we have set it to 
zero. 

The two radii, R±, are both tunable, but with string coupling Qs small or with 
K sufficiently large, R- can be considerably smaller than This is the regime 
of interest for us since i?+ > -R_ typically generate a large hierarchy, be it between 
string scale and inflation scale or between string scale and electroweak scale. 

The local geometry in terms of the physical metric is then 

^ + Rt log{r/R^rf' + dsQ = - ^Rtyf^" {dy' + d4j 

(2.8) 

Up to the change of overall size due to the logarithm, this metric describes a line 
parameterized hj y = — log(r//2+), times T^'^. Together with smooth completion 
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near r = this is called the Klebanov-Strassler throat. The radius of T^'^ varies from 
i?+ near the top of the throat (r ~ to -R_ near the bottom of the throat. At the 
bottom of the throat radius of 5*^ is while radius of 5*^ is shrinking to zero. 

This approximate metric based on the singular conifold must be replaced by the 
one based on deformed conifold as we reach the size of T^'^ to be We will denote 
the value of r there as tq = i?+e~^°, which takes the value 

1 > e'^" = — = e-(^+/'^^-)+l = g-(27rA733,Af)+l (^2. 9) 

For r < To the value of H does not change much, and we then have a redshift factor 
between the top and the bottom of the throat 

1 » H-'/\ro) ^ ^ = |i e-(«V^^-)+i = ^ e-^» (2.10) 

This small factor is responsible for generation of hierarchy between top and bottom 
of the Klebanov-Strassler throat. 

2.2 Hierarchy and the Decay Cascade of Local String Modes 

We identified the cascade of energy transfer from D — D annihilation to SM particles 
as shown in the Figure |21 First part of the process is creation of closed strings loops. 
Then they decay into lighter KK modes. KK modes interact with SM branes to 
excite SM particles. Residual scalar excitations of SM brane decays further into SM 
fermions. Each step of this cascade will be described in Section El Here we give 
qualitative description of excitation modes (particles) in the warp geometry, which 
will be crucial for quantitative estimations of Section |21 

The hierarchy is generated because the total metric has the form 

G = H-^'^g + h (2.11) 

This structure implies that the conserved energy E in the noncompact spacetime 
depends on the position of the quanta in the compact manifold as 

~ H-^'^{y) (2.12) 

This can be seen from the on-shell condition for particle of mass m, 

H^'^g^y^Pu + h'^pipj = -n? (2. 13) 
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where would be the conserved momenta in the spacetime for a Minkowski! metric 
9- 

For degrees of freedom arising from closed strings, the right hand side would be 
represented by the oscillator contribution, and we have 

- 9''''p^.Pu = H-'l\y) (h^'pipj + ^) (2.14) 

with an integer quantized oscillator number N, and the internal Kaluza-Klein mo- 
menta pj. This crude formula is valid only when we can regard the geometry h and 
the warp factor H be both sufficiently slowly varying so that we can regard internal 
geometry h nearly fiat, and H nearly constant. Furthermore in this rough scaling 
argument we also ignore change in string quantization due to fluxes. 

In fact, it is doubtful whether these assumptions are justifiable for most of closed 
string modes we are familiar with in old Calabi-Yau compactification without fiux. In 
the absence of fiux, the Calabi-Yau compactification gives two types of closed string 
excitations. One class is the oscillator modes which we usually ignore for low energy 
purpose since the mass thereof are all fundamental string scale -J^ at least. The other 
class is Kaluza-Klein modes which arise from Fourier analysis of lOD supergravity 
modes on the compact internal manifold. For large volume, these latter modes are 
most relevant. These KK modes are expressed in terms of eigenmodes of various kind 
on internal manifold, which take nontrivial wavefunction throughout the underlying 
Calabi-Yau manifold, and can in no way deemed to be localized in some part of the 
internal manifold except for those with extremely large mass. 

When the fiux compactification involves a warped throat of the kind we discussed 
with large redshift factor, however, a new class of closed string modes emerges. While 
the conifold region is a very small piece of the Calabi-Yau manifold, Q, its associated 
throat geometry is not necessarily small since the physical metric is h = H^^^Q with 
an exponentially large H. Thanks to this, there are both oscillators mode and KK 
modes which are localized deep down the KS throat. Mass scales of such local modes 
are 

rriKK ~ H~'^^^{ro) — , niosciUator ~ H~^^'^{ro)—= (2.15) 



respectively. Provided that somewhat larger than the string scale y/af, it makes 
sense to consider such localized and thus redshifted modes. 
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For these modes at least, the above scahng of on-shell condition is apphcable. 
While curvature of h, gradient of H, and existence of fluxes will change the form of 
the on-shell condition, none of them will have as dramatic effect as the factor 
in front of the right hand side. In terms of the linear coordinate y — — log (r //?+), 
we have 

//-V4 ^ g-s/ (2.16) 

up to a prefactor, and this exponential dependence dominates any kinematics of the 
closed string modes. 

For any massive modes of mass in the unrescaled unit, located deep down the 
throat, we have the on-shell condition, 

-p^p^~e'^°-'W (2.17) 

with the right hand side increasing exponentially as we move up the throat, away 
from conifold toward the bulk of Calabi-Yau manifold. This shows roughly how the 
mass quantization differs between local modes and the rest. For modes deep down 
the KS throat, the mass gap scales either as ~ e~y°/R- or as ~ e~^°/-\/a', while 
everywhere else the mass gap scales as l/R, with R being the hnear size of the 
Calabi-Yau manifold, and 1/Vc/.* 

Therefore, for energy distributed among local modes to escape the throat, one 
must assemble the energy into a few quanta with exponentially large kinetic energy. 
The required kinetic energy must be larger than its rest mass by a factor of » 1. 
Otherwise, the quanta would be simply reflected by a Liouville-like wall. The strength 
of these Liouville-hke walls is dependent on the mass of the particle: the heavier the 
local mode, the stiffer is the wall. Thus, in the presence of this mass-dependent 
potential barrier, localized heavy mode will have tendency to decay near the point i/c 

evo-vc ^ E (2.18) 
m 

to modes with smaller mass scale. 

*Here we are assuming that all moduli are fixed by the flux compactification. For some moduli, 
notably those to be fixed by nonperturbative mechanism, the associated mass scale could be; con- 
siderably lower. However, this separation of two scales associated two types would-be moduli are in 
principle independent of this hierarchy we address, and will be taken to be insignificant. 
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$ = (2.19) 



Any heavy localized mode present deep inside the KS throat would eventually 
decay to lighter modes within the same throat. As the energy cascades down to the 
lighter modes. For substantially larger that a', the lightest modes will be KK 
modes localized at the bottom of the throats, and thus in the intermediate time scale, 
energy will be deposited in these modes. While they are light, they also are massive 
with characteristic mass scale e~y° /R-, and thus are also confined within the throat. 
The associated potential wall is less stiff, and the energy can be stored in a somewhat 
larger volume at the bottom. 

On a quantitative level, what was said above can be seen in terms of eigenmodes 
which obey the oscillator-like equation [201 122| • Here instead of writing the full KK 
mode equation, we rely on a simple massless scalar field eigenmode equation 

Factorizing the eigenmode into outer /inner space parts, and replacing four- dimensional 
Lorentz invariant operator V^V^ by its eigenvalue mj^j^, we get equation for the KK 
eigenmode ^jtikk- This mj^j^ is KK mass as seen by four- dimensional observer, 
'^'kk = — k"^- Taking the approximate form 

H ~ e^^ h = R^ (dy'^ + ds^^ ^) (2.20) 

where i? is a slowly varying function ranging from /2+ to R- which we take to be 
constant effectively. In other words, we can approximate geometry with that of AdS 
warped geometry of Randall- Sundrum type, times the internal compact manifold 
Ti l of a definite size. Singular boundaries of the Randall- Sundrum geometry are 
naturally smoothed out by having the additional internal dimensions. Attachment to 
Calabi-Yau manifold and the cigar-like capping of the bottom, respectively, replace 
UV and IR branes. 

Denoting the quantized and dimensionless angular momentum on Ti i by L^, this 
equation simphfies to 

e'ydye-^ydy + mlj,R\^y - L"] = , (2.21) 

where L? is a contribution from angular momentum. Its spectrum depends on the 
isometry of Ti i. For instance, contribution from 5*"^ sphere will give L"^ = l{l + q — 1) 
where / are integer numbers. 
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As usual with such warped geometry, the equation can be transformed to a Bessel 
equation by taking 

z = e^ ^m^^.L = ^'^™;L (2.22) 

which gives 



z z 



-^m^^^L = (2.23) 



Thus the eigenmodes to this simphfied KK equation are given by hnear combinations 
of Bessel functions 

J^^irnKKRe"^) (2.24) 

with z/^ = 4 + (for L = we shall take combination of functions J2 and Y2). This 
shows that, with the length of the interval in z coordinate of order ~ e^°, the mass 
eigenvalues are quantized in unit of 



e 



-ya 



^rriKK ~ s- (2.25) 



as promised. The mass spectrum is roughly 



e 



-2/0 



rriKK ~ n—— (2.26) 
rC 

with integers n, and for large n, the wave function ^rnKK;L{y) is oscillating near the 
bottom of the throat, while far away from the bottom (small y) is has the asymptotics 

In the long run, the decay will further proceed until energy is mostly retained by 
massless or nearly massless modes of the string theory in question. Assuming com- 
plete stabilization of Calabi-Yau moduli by flux, the only such modes are 4 dimen- 
sional gravity and possibly light open string modes associated with stable D-branes 
which may exist the bottom of the throat. The question is then how the final distri- 
bution of energy will look hke between the bulk gravitational sector and the localized 
open string sector of stable D-branes. 

One important point in pursuing this question is that the only light degrees of 
freedom is the 4 dimensional gravity, but this couples to the modes localized at the 
bottom of KS throat very weakly. The redshift factor H~^/^ reduces the effective scale 
of energy-momentum by an exponential factor and pushes down inflation scale and 
subsequent reheating scale as well. This is essentially the physics of Randall- Sundrum 
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scenario, realized in string theory setting, and can be understood from the fact that 4 
dimensional graviton has the wavefunction profile H~^/'^{ii) in the internal direction. 
Any localized mode at the bottom of KS throat will have very small wavefunction 
overlap with 4D graviton and thus cannot generate much gravitational energy. 

On other hand, open strings on D-branes will couple to KK modes without this 
exponential suppression but there still is a volume suppression if they are for example 
on D3 branes transverse to the compact Calabi-Yau. In order to see how much 
energy is deposited to what species of particles, we must pay more close attention to 
interaction at the bottom of KS throat, which will be discussed in a later section. 

2.3 Subtleties 

The above equation ()2.23j) is usually obtained via vast simplification of actual KK 
mode equations. In particular, the capping of the throat near y = i/q is not faithfully 
reflected, which leaves determination of precise spectrum difficult. On the other hand, 
the robust part of above estimate is that the lowest energy is of order ~ e~^°/i?_ and 
also that subsequent gap between adjacent eigenmodes is also ~ e"^"//?-. 

An important subtlety to bear in mind is how far one should trust this equal- 
spaced spectrum. This linear analysis suggests that a throat of length R^yo has an 
exponentially large number of states of order 6'^^° due to its very small massgap. If we 
take some state of mid-range value of the mass m' such that 1/-R+ ^ m' ^ e~^°/-R_, 
number of states it can decay to is of order (m'i?_)e^° and this will induce a very 
large width to the eigenmode thus obtained. For larger enough m', it is therefore 
reasonable to expect that this linear analysis is misleading. 

This problem also manifests in the shape of the eigenmode. While it takes a simple 
innocuous form in the conformal coordinate z, its behavior in physical coordinate y 
is much more drastic for high- lying modes. With mR-e~'^° ^ 1, the Bessel functions 
oscillate wildly near y = yo, and its derivative could be larger than the local string 
scale. So once we have ttixk larger than the local string scale, it is unclear whether 
supergravity approximation can be trusted. 

Thankfully, however, analysis of the present work is not be affected by this am- 
biguity. We will consider closed string oscillator modes with perhaps up to 100 
oscillators, which eventually will decay to low lying KK modes. As we will see later. 
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energy deposited to KK modes will quickly thermalize among themselves, and since 
the energy scale of initial state right after inflation is or order ~ l/An^gsa'"^ < 
relevant KK modes are those below string scales. We never rely on very high scale 
KK modes who precise nature would need more careful analysis. For bulk estimate 
of reheating processes, even details of low lying KK modes does not enter other than 
their numbers. 

One important exception to this is how the lowest energy eigenvalues depend on 
the angular momentum L^. Later we will consider as approximately conserved 
quantum number in the full Calabi-Yau compactification, this precise spectrum at 
low end could be important for lifetime estimate of long-lived KK mode carrying such 
quantum numbers. We leave it to future study. 

3 Decay of D-Branes at the End of a Brane Infla- 
tion 

One attractive class stringy inflation models involve unstable D-brane system 
whose elevated vacuum energy drives the inflation While introduction of un- 

stable brane system is an novel element, this difference does not seem to generate 
much new flavor in terms of studying 4 dimensional low energy effective theory during 
inflation. As is typical with inflation, maintaining sufficient number of e-folding and 
at the same time having a graceful exit is not a small problem, and these are just 
translated to more geometrically constraints on the underlying string theory. 

When it comes to reheating process after the end of inflationary era, however, 
behavior of brane inflation could be very different from ordinary field theory models. 
Sometimes the so-called tachyon effective action is invoked, but one cannot take 
this tachyon effective action too literally. While this low energy approach has been 
immensely successful, mathematical results one find of it must be reinterpreted with 
care. For instance, the so-called tachyon matter is known to survive the decay process 
and takes up all initial energy in the unstable brane system, and behaves like a perfect 
fluid of very massive noninteracting particles. However, the system in question started 
with an open string description which should not be valid by the time D-brane has 
decayed. This so-called tachyon matter turns out to be a coarse-grained view on 
the underlying physical state, namely a certain distribution of highly excited closed 
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string states. 

In the conventional inflationary models, reheating question centers on how effec- 
tive an inflaton decay can excite other degrees of freedom. Here, the situation is 
reverse. The initial reheating process is such that 100% of the energy density respon- 
sible for the inflationary phase is converted to heavy degrees of freedom that have 
nothing to do what drove the inflation. The right question to ask here is how this 
huge amount of energy density is eventually distributed among different hght degrees 
of freedom. Since the initial stage of decay produces a lot of heavy closed strings, 
there is an inherent danger of closed strings dominating the process. Assuming some 
kind of brane worlds scenario for standard model sector, this would be a disaster for 
brane inflation to be viable, or looked backward a very efficient and simple tool for 
eliminating many stringy inflation scenarios. 

In this section, we will describe a very efficient and viable reheating process for 
the case of single throat scenario. In later section we will discuss under what circum- 
stances multi-throat scenario can offer a viable reheating process. 

3.1 D — D Annihilation and Closed Strings Production 

The end point of inflation is D?) — Z)3 pair annihilation, the step 1 at the Figure El 
Complete description of non-BPS system in string theory is a complicated problem 
|28j . The hallmark of the D-brane decay, as opposed to decay of inflaton in fleld 
theoretical models of inflation, is the fast and complete conversion of unstable D- 
brane energy into massive closed string modes EUl E]- This is what we flnd 
from tachyon dynamics in low energy [H2l IHHl IHH IHH] but also supported by stringy 
computation |2ZlEHlEni using decaying boundary state jlHl EH 112 ESI • 

The pair of isolated — Z)3 brains annihilates into excitations of close strings 
loops with the average energy in the bosonic string theory [38^ E is 

^^Y. I d^'^'^'^k^ D{N) e-^^-iv.fe ^ (g^) 

where ujn,k = \/k'j_ + 4A^, V3 is the three dimensional volume of the branes, k± is 
the momentum of the closed strings transverse to the branes, D{n) is the number of 
closed strings oscillator states, for large 
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Cascading Energy from Inflaton to Radiation 




Figure 2: Identifying the channels of D-brane decay 
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in bosonic theory q — 27/4, d — 25; in superstring theory d — 9 and q is not known. 

It is interesting to compare this result with the inflaton decay in the QFT tachy- 
onic preheating, which is characterized by the high occupation numbers. In brane-anti 
brane annihilation, the product states are characterized by two numbers, oscillator 
state n and momentum k. analog of the occupation number would be e~^'^'^^''= which 
is less than unity. However, QFT has no analog of the oscillator number. If we sum 

— * 

over n for a given k, in principle we can get the number exceeding unity. 

With a view toward reheating from brane inflation, we will not need much of the 
details of the decay. Let us summarize the main characteristics of closed strings from 
decay of unstable D-branes: 

Initially massive oscillator modes are produced. The probability of a particular 
closed string mode with energy u;N,k to be produced scales as e~^'^'^^>'=, but this expo- 
nential suppression is exactly cancelled by the Hagedorn growth of density of state 
D{N) at larger energy. The upshot is that for each oscillator level N, roughly the 
same amount of energy is deposited. 

Since boundary state is formulated at — > limit, one must introduce cut-off 
to emulate backreaction of the boundary state to production of closed strings. For 
unstable DO, this natural cut-off is rus/gs with the local string scale m^, and with 
this cut-off the energy in the produced heavy closed string account for all energy in 
the unstable D-branes. The highest oscillator stare which is expected to be excited 
is estimated from uj^ax ~ l/s's- For the string coupling Qs ~ 0.1 we get A^rnax ~ 100. 

The probabihty distribution e~^'^'^^>'= also implies a narrow velocity dispersion. 
The usual string on-shell condition uj ~ ^Jp^ -\- ml N, implies that the transverse 
velocity of typical closed string mode is at most of order ~ -y/rn^/ y/m ~ ^/^. Thus, 
kinetic energy of the produced closed strings tends to be smaller than its rest mass. 
The string loops are non-relativistic, only very slowly can move around their birth 
place in the volume Vds- 

If the unstable D-brane system decays well inside the KS throat with a redshift 
factor if~^/^(ro) and a large radius all of above should remain true quahtatively. 
The main difference is in the string scale m^. Since we are discussing energy in terms 
of 4-dimensional metric g, the string scale that appear above should be related to 
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the fundamental string tension by 

„j.£;;^!M.,-.^ (3.3) 

Thus, the decay of D-brane at the bottom of a KS throat produces strings of mass 
m and energy E of order e~''^^ / {qs^/o') at most. 

Prom the four dimensional perspective, energy density of the closed loops is the 
same as the energy density of the £)3 — 53 pair. We have to take into account redshift 
of the brane tension in the warped geometry by factor e"^^" 

which is the scale of the warped inflation. 

The timing of — D?> annihilation we estimate as 

Ail ~ e^°v^ . (3.5) 

Strictly speaking, the string computation was done for decay process involving 1- 
point function of closed strings, and for this reason can be rehed on only for unstable 
DO or unstable D-branes wrapped on a small torus. The same computation for more 
generic cases is inconclusive because the 1-point emission is not the dominant decay 
channel. In fact, energy output from such 1-point emission can be computed and 
can be shown to be well short of the expected energy output. Mathematically this 
happens because the produced string mode cannot carry momentum longitudinal to 
the D-brane. However, this is an artifact of 1-point emission processes. 

For higher dimensional cases, the dominant process should involve simultaneous 
production of two or more closed strings, which has more phase space volume along 
the longitudinal direction. Although such a multi-point amplitude would be sup- 
pressed by powers of string couplings, the gain from phase space volume can easily 
overcome this. For instance, consider an unstable Dl which is infinitely extended. 
Two point decay would be suppressed by one more power of Qs, but the available 
phase volume grows linearly with the energy of the created pair of particles. Intro- 
ducing a natural cut-off I/qs in string scale, the phase volume will be of order I/qs- 
This cancels the coupling suppression from the string diagram easily; we expect this 
process to dominate the decay process with the quantitatively same characteristics 
of produced closed strings. 
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Consideration of low energy approach using tachyon effective action has been 
developed independent of the string theory computation and tells us pretty much 
the same story. The so-called tachyon matter [lUJ HI] and string fluid [321 ESI EH 
El] , which emerges from this low energy approach, have been studied in depth and 
compared to closed string side |8T| I45j . for unstable D-branes of arbitrary dimensions. 
The result shows that all energy is converted into heavy closed string excitations, 
possibly together with long fundamental strings. 

3.2 Closed Strings Decay to Local KK Modes 

Now that we identified the initial state right after the end of brane inflation, we 
must consider the subsequent decay of massive closed string thus produced. This 
corresponds to the step 2 at the Figure El 

The local closed string modes will then decay to whichever are lighter degrees of 
freedom around. There are two types of light modes at the bottom of throats. One 
are KK modes, which would be lighter that oscillator modes, as long as /?_ is not 
too near the string scale. Another are light open strings modes associated with extra 
D-branes. One easy way to achieve this is to have some extra D3-branes surviving at 
the end of annihilation. While decay amplitude to a pair of open string modes would 
carry one less factors of g^"^, it is also suppressed by volume effect ~ (a'/i??.)^/^. 
Even with R^j\foi' slightly larger 1, the latter effect will easily compensate for the 
former, and will favor decay to KK modes first. In this work, we assume that 
is large enough to justify field theory and gravity analysis employed here. Even if it 
is necessary to extrapolate to smaller value -R_, our guess is that basic qualitative 
estimates we offer in this paper would remain valid. After all, what really matters at 
the end of day is that energy settles in some open string sector that would contain 
the standard model, and non-SM dangerous relics. 

In the fiat background massive closed string with oscillation number n decays 
into two closed strings with oscillation numbers N — N' and A^', A^' < A^, with the 
coupling strength gg. Subsequently, N — N' and A^' states further cascade into states 
with the smaller oscillator numbers. The final states will be KK modes, fraction of 
which is in massless gravitational radiation in the bulk. Gravi-tensor projection of 
the bulk gravitons to four dimensions give describe usual gravitational waves in four 
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dimensions. Therefore, in the inflationary scenarios which have D — D annihilation 
at the end, we shall address the problem of overproduction of gravitational radiation. 
One of the possibility to overcome the problem is to arrange annihilation not between 
a D—D pair, but a system of D and several D branes, to dilute to amount of gravitons. 

However, decay of closed strings in the warped geometry is significantly different 
from that in the fiat background. Although the energy of the unstable D-brane system 
is all deposited to the highly excited strings, none of these modes can overcome the 
potential barrier toward the top of KS throat, since there is factor difference 
between strings scales between the top and the bottom. 

The only exception to this would be the part of energy radiated into massless 
modes. In the current setup the only massless graviton is the 4-dimensional variety, 
whose wavefunction comes with the same exponential suppression at the bottom of the 
throat. Another way to say this is that the effective tension of the unstable D-brane is 
redshifted by a power of 6"^° relative to the Planck scale which weakens its couphng 
to the 4 dimensional gravity. It also can be seen in the following way. The wave 
function of massless four dimensional graviton (zero KK mode), after rescaling with 
respect to the warp factor, is homogeneous relatively to inner dimensions. Therefore 
decay of the closed strings into gravitons is suppressed by the four dimensional Planck 
mass, while decay into massive KK modes is suppressed by the redshifted string scale. 
Therefore, we have very specific effect of the warped geometry of the exponentially 
suppressed production of gravitational waves from D — D annihilation. To the first 
approximation, no energy is deposited to the 4 dimensional graviton directly. More 
specifically, the energy fraction deposit in gravitational waves radiation is e~^^° of the 
total energy of radiation. This figure depends on the energy scale of infiation. Recall 
that amount of energy density in gravitational waves at the moment of BBN shall be 
no more than several percent of the total radiation energy density. 

Apart from the 4 dimensional gravitons, the coupling to which is universally small 
deep inside the throat, the next fight degrees of freedom are local KK modes. To see 
this clearly and also for a later purpose of computing branching ratio, let us expand 
the 10 dimensional Einstein action according to the compactification. Starting with 
10 dimensional action. 




(3.6) 
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with = (27r)''Q;"^/2, take as before G — H~^/'^g + h but allow small fluctuations in 
g and h. The former gives the 4-dimensional gravitons while the latter is prototype 
for KK modes from the compact manifold. Expanding the 10 dimensional action and 
keeping terms quadratic to the small fluctuations, the first term 



IGttG 



/ dx'-^'V^R, (3.7) 



1 ^ V(i 



is the Einstein action in 4 dimensions with the effective Newton's constant Gn such 
that 

j dy^VhH-^/^ (3.8) 

"10 

Note that the KS throat region contribution to ~ Vq is of order R^, so we have a 
bound Vq > R^. Large i?± will necessarily make the volume of the compact manifold 
somewhat larger than fundamental string scale, and weaken gravity further, regardless 
of the exponential redshift factor. 

In our convention, the gravitational redshift due to the warp factor is manifest on 
the matter effective action. If we expand h = hxs + the remaining terms has the 
following general structure, 

[ dx^+^yf^d dy''y/^H-^'''5hVl5h+ I dy''VhH-HhS/l^Jh] (3.9) 

which are the kinetic term and the mass term, respectively, of internal KK modes. 
Further expanding 5hva. terms of eigenmodes localized at the bottom of the throat, 

5/^ = E'^&(^MV) (?/) + ••• (3-10) 

where 5h{n) denote local KK modes, all supported at the bottom r ~ tq. 

As a matter of convenience we will normalize 5/i(„) to produce canonical kinetic 
term for 0^']^, so that these second class of terms produce, 

\ f dx^+^v^E ((V4t)^ - {milm^K?) (3.11) 
The KK masses m^]^ of the locahzed modes are of order 

mi^^e-yi-^^ (3.12) 
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where we introduced y{n) whereby Sh(^n) occupy region at the bottom of the throat 

y > y{n)- 

As we saw above, the initial product of D-brane decay is numerous heavy closed 
strings of mass ~ e''^°/gsVc^ or smaller. With not too small local KK modes are 
much lighter and more numerous than these local oscillator modes, and will couple 
to the closed string modes in usual 3-point diagram to these oscillator modes with 
the coupling being essentially Qg. The KK modes thus produced can have energy or 
mass up to ~ e~y° /qsVc^, which implies that the energy is distributed in KK quanta 
and excites the throat up to yo — 2/ ~ ^og{R-/ gs\fc^)- 

In contrast, the 4 dimensional graviton, which is the only light mode that are not 
localized at the bottom of the throat, has much higher mass scale of its own, 

_ _ 1^ I Rl e f27nMKY^' e 



which is exponentially larger by a factor of than that of other scattering processes 
among local string modes and local KK modes. Here the numerical constant e which 
came from the definition of k\q is 

/ 9 \ 

(3.14) 



Up to this stage, the primary decay channel of closed string produced from the D- 
brane decay will be into local KK modes. 

Finally we estimate the timing of closed strings cascading into light KK modes. 
Rate of decay of the individual close string is about F ~ g'^e~^°/Va'- The longest 
process is the closed string with the oscillator number N^ax goes through N^^x decays, 
so the upper bound on the net decay time is t^t^ ~ ^maxl^- Since N^^ax ~ l/fl's; 
estimate 

Ai2 ~ i/7V° ^ . (3.15) 

his timing is longer than the brane annihilation KK mods are ultrarelativistic. 
Indeed, closed strings have energy rriciosed ~ ^ J^ , while KK modes have the mass 

rriKK ~ # so that ^^^^ ~ 4=7 -• 

Next we proceed with the decay of KK mode into open strings modes, which can 
be associated with the SM particles. 
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3.3 KK Modes Decay to Open String Modes 

Let us note that upon exciting KK modes, the corresponding metric perturbation 
along the compact direction has the size 

gsKwH^/^iro) („) 
^Ky) ^3 ^KK (3-16) 

for each KK eigenmode with the canonically normahzed 4 dimensional massive fields 
4>iiK- This is what couples to open string degrees of freedom directly as we will see 
shortly. With this, let us consider how KK modes couple and decay to open string 
modes on a D-brane transverse to the Calabi-Yau direction. 

One may be alarmed to see the exponentially large factor if^/^ in this expression, 
but all this does is to rescale the dimensionful parameters in terms of the local 
redshifted scale, since 



The mass scale of relevant (pxK^ are anywhere between e~^°/i?_ and e~^"/v/a', so 
as long as R^ is not too small in string unit, the fluctuation is small and treatable as 
a perturbation. 

With this in mind, let us consider a probe D3 brane located at the bottom of the 
tip. The Born-Infeld action can be written as, 

- g^3^ / dx^+^H-\r)^-V>ei {g^, + 2T:a'F^, + Hy\r)hijd^Yi d^Y^) (3.18) 

where r = t{Y) and Y^{x) represent transverse fluctuation of D3 brane along the 
compact directions. The primary interaction between 5h and D-brane appears from 
the leading expansion, 

^J^^ j dx^^'^g {h-^/\Y) {^-hijd^Y'd^Y'g^^^ - H-\Y)^ (3.19) 

Here we have chosen the longitudinal coordinate system such that 

Qnu = r]^^ + 5g^,^ (3.20) 

where 5g represents the 4-dimensional graviton. For the moment we set = so 
that g = T). 



23 



In addition, there is another potential term from the minimal coupling to the 
background R-R 4-form potential C4. For this let us recall that KS geometry comes 
with 3 types of fluxes; NS-NS 3-form flux — dB^, R-R 3- form flux F3, and 5-form 
flux which is related to the previous two by 

F5 = dCi + B2AF3 (3.21) 

and is constrained to be self-dual. A minimal C4 for this self-duality constraint to 
hold can be determined from B2 and F3. Relying on the explicit solution of Ref.[KS], 
we find 

C4 = cH-^ dx^ A dx^ A dx^ A dx^ (3.22) 

D3 branes have the usual minimal coupling to C4. The proportionality constant 
c is such that this contribution from R-R coupling cancels the tension term from 
the Born-Infeld piece exactly for a probe D3 brane and double it for probe anti-D3 
brane.^ Thus, the potential term would be either cancelled or doubled in actual 
dynamics. 

Expanded up to quadratic order in y, and canonically normalizing Y fields to Y , 
we have the following general form of action 

1 / y=i ((fiy )^ - ^^y ^ + May?) (3.23) 

where we ignored a multiplicative constant of order one in front of the last, cubic 
interaction term. The mass term is /i/)3 ~ e^^''/i?_ for anti-D3 brane, but for D3 
brane, /i^s = within the present framework. 

With a single throat containing both inflation and standard model, there should 
be a further hierarchy generating mechanism. One such possibility is to have super- 
symmetry unbroken to much lower energy scale. For this, we must then consider 
leftover D3's rather than anti-D3's, for otherwise we would end up breaking super- 
symmetry at the inflation scale. The story of SM brane phenomenology is out of the 
scope of this paper, we just assume small but non-zero fi which is larger than TeV 
scale. 

^This is the same cancellation/doubling that is usually employed in deriving a slow role potential 
for a D3-anti-D3 pair. 
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In fact, it is not even necessary to assume that the transverse scalar fields are 
the main decay channel for the KK modes into SM brane world. The above three- 
point coupling is generic enough to work on any low lying degrees of freedom on the 
SM brane world. For instance, worldvolume fermions would also have such a 3-point 
coupling, and can absorb energy from KK modes directly. Since we are only interested 
in rough estimate of reheating process, we will work with the scalar Y and its specific 
form of coupling while keeping in mind that the result applies to all standard model 
sector. The interaction between (pxK and Y has the general form 

1 



A 

with 



(I^KKidYy (3.24) 



1 f7l9',M'^'/' 



e^'o — = ( — ) e^o Va' . (3.25) 



A Ri \ 32 ^ 

This type of interaction is typical for interaction of the radion in the brane-worlds 
scenarios [46 . The rate of decay of KK modes into Y is 

For anti-D3 brane we have a problem of decaying into Y particles, since /i£)3 and 
niKK are comparable. 

On the other hand the mass /i^s can be arranged to be smaller than tukk, and 
therefore D3 branes as the SM brane is preferable. We proceed assuming Thus the 
timing of decay of KK modes is At-^ = I/Tkk 

Atg = 327r {^] = 327r ^ (3.27) 



^ 32 ^ 

If we adopt R- to be larger than y/a', then so far this is the longest process in the 
chain of energy transfer from inflaton field. As we encounter many times throughout 
the paper, once again our conclusion depends on specific range of parameters. Our 
choice here is based on justification of computation in low curvature KS geometry. 
As we mentioned above, KK modes can also decay directly into SM fermions, with 
the rate of decay similar to (j3.26p . 

Before proceeding with the decay open string modes into SM particle, we shall 
discuss separately a special case of KK modes, which can be long living and most 
dangerous for the whole scenario. 
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4 Thermalization and Dangerous KK Relics 



In the previous section we consider decay of massive (but light) KK modes into 
excitations of the SM brane, and it was found that in the model with radius of KS 
throat around its tip i?_ significantly larger than ^/d (as expected in the supergravity 
description of KS geometry), decay time of KK modes (j3.27|) is significant. 

KK modes are self-interacting. Therefore let us check if KK modes are thermal- 
ized or not before they decay into open string modes. In this section we address 
thermalization of KK modes. Then we consider residual decay of open string modes 
into SM particles, and thermalization of SM particles. Most importantly, we iden- 
tify the problem of long-living KK modes, which is a serious problem for the string 
inflation scenarios. 

4.1 Thermalization of KK modes 

To check whether or nor KK modes are thermalized we have to compare the time 
of relaxation of KK modes towards their thermal equilibrium, r, and time of their 
decay Ata. 

Relaxation time towards thermal equilibrium is 



where n is the 3d number density of KK modes, a is cross-section of their rescattering 
and V is their typical velocity. We have to estimate each factor in ()4.1|) . 

In principle each factor in (j4.ip depends on the expansion of the universe. We 
however argue it can be ignored. Indeed, let us estimate the value of the Hubble 
parameter H immediately after brane-anti brane annihilation. From the Einstein 
equation 



The inverse Hubble parameter, which is a typical time of expansion, is suppressed 
by the small factor e"^^'', while all time intervals of interaction including r, as we 



1 



(4.1) 
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will see below, are suppressed by factor e~^" only. Therefore here expansion can be 
ignored. 

Number density can be estimated by ratio of the total energy density e and energy 
per KK mode, e can be taken to be the energy density after the brane annihilation 
(13. 4j) . and energy per KK mode to be comparable with the energy of closed string 
loops e'^^'/v^, from which they originated. Thus 

g-3yo 

" ^ 4^^W^ ■ (4-4) 

KK modes have masses ttikk ~ while they are created from closed strings of the 
mass %=■, so that KK modes are relativistic v = 1. 

More elaborated is estimation of a. Interaction of KK modes can be derived in 
the following way. Again, take lOd action ()3.6|) with the ansatz G = H~^^'^g + h as 
before in Section The four dimensional action, extended up to non-linear terms 
with respect to h, contains the term -gr / d^^^x^J—g{hVhY. After decomposition 



h = hxs + it contains, in principle, three legs and four legs interactions, which are 
comparable. We take four- leg interaction {6h'VShy. Cross section of this interaction 
a is the ratio of squares of coupling and energy, the coupling is ^"'^^"2°^ , thus coupling 
is rather strong ~ g'^. This is because KK modes are very energetic, k ~ ^y=-, so 
that momentum factors compensate the suppression by the local string mass. The 
estimate of relaxation time is 

47r3 



-e^oVa' • (4.5) 
9s 

Relaxation time of KK modes rescattering is much shorter than the time of decay of 

unstable KK modes At^. Therefore we can treat all KK modes as particles in thermal 

equilibrium. We encounter rather unusual situation in cosmology when 

non-SM particles are set in the thermal equilibrium first, even before SM particles are 

produced! 

Therefore the issue of the reheat temperature is split into two issues: what is 
the reheat temperature of KK particles, and what is the reheat temperature of SM 
particles. 

To calculate reheat temperature of KK particles, T^k, we have to convert the 
energy density KK particles into thermal energy of KK plasma. Energy density of 
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KK particles 

e~c^, (4.6) 

where numerical coefficient c depends on equation of state, for radiation equation of 
state of KK particles c = |. The end of inflation fixes the initial value of t = tg 
in ()4.6|) . We have to compare ()4.6|) energy density ()4.6|) at the moment t = to + 
with thermal energy of KK particles. Suppose KK particles are thermalized being 
relativistic. Then their energy is 

e = i^9kkTkk , (4-7) 

where qkk is a number of degrees of freedom of KK particles. 

A subtle point is that to ^ t. Indeed, the end of infiation to is defined by 
equalizing ()4.6|1 with the energy- density of original D3 — D3 pair ()3.4|1 

or similarly from to ~ ^/H with H from ()4.3p . We find 

to ~ J 4:71^ eg, e^y" ^ (4.9) 



1/Va' 

thus to ^ r. 

Therefore, comparing c-r^ with (14.71) . we have 

g-s/0 X 



1/4 1/4 

Thus, reheat temperature of KK modes up to a factor gj^j^ is comparable with the 
inflaton energy scale. 

It is instructive to compare Tkk with the mass of KK particles rriKK ~ 
Unless gKK is very large, reheat temperature of KK particles is larger than their 
mass, so they are thermalized as relativistic particles. A new interesting element here 
is gxK factor. As we seen in Section 3.2, KK modes have energies up to e^y° / g, ^fo' ■, 
the number of degrees of KK modes gxK of the mass rnxK is then proportional to 
R-l ga\fa' . However, formula ()4.7|) is valid only for relativistic KK modes, i.e. for 
modes with masses less than T^K. For non-relativistic KK modes we have to use 
formula exK = ruMM .tikk with tikk = (^^^ff^)^''^ g"""^^/'^. 
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4.2 Disappearance of Open String Relic 

Among the product of decaying KK modes into open string sector of SM brane, 
there was the scalar Y, describing displacement of SM brane. Specific field content 
of such scalars are dependent on how SM is realized, but nevertheless we must worry 
about such particles as undesirable moduli fields in cosmological sense. We need 
the mechanism to give rid off Y, even if such a field exists on SM brane world. 
Fortunately, the gauge theory associated with the brane contains fermions ip in the 
adjoint representation. Decay of (pKK can go also into fermions, we are interested here 
in the troublemakers Y. The scalars Y interact with the fermions with the three- 
linear coupling ^J^gipYip. Therefore Y shall completely decay into lighter fermions, 
the rate of decay is 

r(F-^V^) = ^. (4.11) 

Remember that the mass ^£,^ here is introduced on the phenomenological ground 
beyond the brane inflation scenario. For /i£)3 exceeding Tev scale cosmological con- 
straints are satisfied with a large margin. Thus decay of the scalars Y by itself is not 
essential process of the reheating, as far as SM phenomenology from string theory is 
successful. 

4.3 Reheating of SM Particles 

Let us calculate the reheat temperature of the SM universe in the scenario. The 
longest process in the chain of energy transfer is the decay of KK modes into open 
string modes given by time At^ ()3.27|) . Now we can compute the reheat temperature 
of SM sector after the warped brane inflation. 

For this we have to convert the energy density stored in the thermal radiation of 
KK modes 

. ^JMk^ (4,2) 



into thermal energy of SM plasma 



4(to + Ahf 



esM = ^9*T^ . (4.13) 

Again, as in the Section 4.1, we shall compare to from (j4.9|) with the decay time 
of KK modes At^ from ()3.27|) . The answer is model- dependent. In principle, i?_ 
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can be as small as string scale, see e.g discussion in ^3]. In this case At-^ can be 
smaller than and estimation of SM particles temperature is similar to that of KK 
modes of Section 4.1, just with replacement of Qkk- However, in the toy model of 
KS geometry in supergravity limit with i?_ > 10^^, for inflation at energy scale 
lO^^Gef with e~^° ~ 10^^, we have At^ greater than to by a large margin. 

Energy transfer completed when the Hubble parameter H{t) ~ 1/t drops below 
the rate T^k, i-e. when t in (j4.13|) is equal to At^. In (j4.13|) g^: is the number of SM 
degrees of freedom, ~ 100. We have 



This estimation is not sensitive to the equation of state for KK modes, if they would 
have the matter equation of state, it changes only numerical coefficients in (j4.12p . 
Reheat temperature ()4.14p is suppressed by two small factors, the redshift factor 
e"^"/^, and by the ratio of scales Reheat temperature of SM sector is lower the 
reheat temperature ()4.10j) of KK particles 

Suppose the energy scale of brane inflation at lO^^^GeV^, e^^" ~ 10"*^, and geometry 
is such that i?„ ~ lOv^- We get ~ 10''G'ef . Now suppose the brane inflation 
is the low energy Tev scale inflation with e~^° ~ 10^^^, and ~ lOv^. Then 
Tr ~ lOOGef , which is cosmologically acceptable. If we choose however, i?„ < 
lOO-s/cy, reheat temperature drops below Mev scale, the temperature of BBN, which 
is unacceptable. Reheat temperature in the scenario is generally lower than that in 
QFT inflation. Therefore we conclude that 

the space of warped brane inflation is constrained from the lowest value of reheat 
temperature. 

4.4 Problem of Long-Living KK modes 

Above we silently assumed that KK modes all decay into open string modes through 
three-legs interaction ()3.24j) which make the decay complete. This however is not 
true for a special class of KK modes for which three-legs interactions like ()3.24j) 
are forbidden. Now we are coming to the critical problem of the whole scenario we 
identified so far, the problem of long-living KK modes. They are massive so that after 




(4.14) 
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universe expansion dilutes energy density of radiation they become dominant. Unless 
extreme fine-tuning, these particles will have unacceptably large energy density flxK- 
In the particular model of warped inflation with KS geometry, there is exact 
isometrics and S'^ of six dimensional interior manifold at the bottom of the 
throat. Suppose more generally we have S*^ isometry inside X^, and 6a are angular 
coordinates of 5*"^. KK modes will be described by the harmonic expansion of rank 
two tensor with respect to the eigenmodes SH'^^^Oa) 

VaV^/i^^) = /i(^) , (4.15) 

(we drop 5 from 5h for transparency), with the Laplace operator on 5*^ in the right 
hand side of the equation. The eigenmodes 5h'^^^ are characterized by the set of 
conserved quantum numbers, associated with the conservation of angular momentum 
and its projections. For example, for the eigenmode of spin s and with isometry 
^0(4) ^ SU{2) X SU{2), we have EH 
L2 = (/ + 1)2_(s + 1). 

The amplitudes of interactions of KK modes include three-legs decay interactions 
with the factor / dVtq h^^^ which vanishes. The only non- vanished amplitudes are the 
annihilation and inverse annihilation amplitudes / d^lqh^^^h^^'^ which conserve the 
angular momentum. 

In the KS geometry, around the bottom of the throat we have isometry of 
sphere with the radius As far as we place SM brane at 5''^, the whole 5*^ isometry 
is distorted by the presence of D3 brane, and some quantum number of the original 
5*0(3) isometry are not precisely conserved. However, residual 5*0(2) isometry is 
intact. Another two directions have isometry of 5^ but its radius shrinks. This 
would blow up corresponding KK mass, so that corresponding KK modes will stay 
away from the bottom. 

Rescattering of KK modes with angular momentum brings them into the thermal 
equilibrium with the rest of KK sector. However, in an expanding universe the number 
density of KK modes with angular momentum freezes out. Indeed, massive KK 
particles sooner or later becomes non-relativistic. Number density of non-relativistic 
particles nxK = ( "^'2^'^ ) e""^^^/^ is decreasing exponentially as temperature T 
diluted with expansion, rate of annihilation and inverse annihilation is exponentially 
decreasing, and when it becomes comparable with the expansion rate of the universe. 
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the abundance of long-living KK modes freezes out. 

The problem of KK modes with angular momentum in KS geometry has similarity 
with the problem of heavy KK modes in supergravity noted in |24j- 

Suppose the angular KK modes are stable. Then the figure of the primary interest 
would be the present day abundance Qrk of the stable KK modes. To estimate Qkk, 
we have to use the standard theory of freeze-out species in expanding universe, see 
e.g. llZj. In principle, particles may freeze out being relativistic or non-relativistic, 
depending on the details of the model. Simplest estimation of Qkk is for relativistic 
freeze out. In this case we have 

Qkk ^0.16 — (4.16) 

gKK eV 

where gst and qkk are the numbers of degrees of freedom of freezed out stable KK 
modes and all KK modes correspondingly. The value of Qkk is of order 10^^ for the 
10^^ GeV scale inflationary throat. 

However, the isometries one find in KS throat is only approximate in the context of 
Calabi-Yau compactification, and therefore there is no absolutely stable KK angular 
modes.. While a KS throat that extends infinitely to UV region do have exact SU{2) x 
SU{2) isometry, the actual internal manifold involves cutting-off the UV region by a 
finite Calabi-Yau manifold. This has the effect of distorting the small r part of the KS 
throat and destroy the isometries. This then propagates toward IR end of the throats 
in such a way that, even at the bottom of the throat, the angular momentum quantum 
numbers are not strictly conserved. In dual field theory language, the attachment to a 
compact Calabi-Yau induces ceratin global symmetry breaking perturbation. Ref. |1H] 
estimated the leading supersymmetry preserving operators to be of dimension 7. 

This then translates to the typical width of the lowest angular momentum mode 
to be 

m{m/Mpf (4.17) 
or equivalently to lifetime of such long-living relic 

■MpV 1 



m J Mp ^ ^ 

With an inflation scale m lower than 10^^ GeV this could easily cause a problem, 
since KK modes live longer than 100 sec, their energy release destroys BBN. If they 
live longer, then there is too much dark matter content. 
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Furthermore, this problem is potentially much more severe if the standard model 
is realized in another, longer throat with such approximate isometrics. As will be 
described in next section, energy transfer from the inflation throat to a standard 
model throat occurs via quantum processes of tunnelling and oscillation, the other 
throat with standard model is likely to be of much lower energy scale. The energy 
transfer then will involve highly excited KK modes in the standard model throat and 
thus will produce relatively large amount of the such angular KK modes. With longer 
throat, the approximate isometry is more protected, and lifetime of such angular KK 
Modes would be significantly longer. 

We should warn the readers that, most likely, existence of such approximate isome- 
trics is not generic. It is true that all known Sasaki-Einstein manifolds which could 
play the role of Ti^i of KS throats are equipped with some angular isometrics, if 
smaller than that of Ti^i. But this is probably result of limited technology on our 
part in finding explicit examples. In dual field theory language, such global symme- 
tries are not required supersymmetry in any intrinsic way. It is easy to envision that 
generic throats in type IIB flux compactification involves no such approximate isome- 
trics at all. Nevertheless our present consideration should exclude certain subset of 
such inflationary scenarios involving relatively symmetry throats. In particular, the 
face value model based on KS throat has problem of long-living angular KK modes. 

5 Issues with Multi-Throat Scenarios 

So far we studied in some detail how remnant of D-brane annihilation cascade down 
to standard model sector, assuming that a standard model is realized at the bot- 
tom of the inflation throat as a brane world. Apart from potential problems with 
long-lived relics associated with approximate isometrics at the bottom of the throat, 
realistic reheating seems possible, although its detail differ from conventional reheat- 
ing process. For one thing, the reheating process involves two distinct phases, where 
localized string modes and KK modes are first created and thermalized and then later 
this energy is transferred to open string sector and thermalize at lower temperature. 

SM throat requires the choice of warping factor e~^° ~ 10~^^. On the other 
hand KKLMMT model suggest inflationary throat to have e~^° ~ 10~^ to have right 
amphtude of cosmological fluctuations. This is one of the motivation to consider the 
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multiple throat scenarios. This constrain can be relaxed if we admit another source 
of cosmological inhomogeneities related to the modulated cosmological fluctuations 
|36j . Notice also that the scale of inflation is tightly related to the SUSY breaking 
scale, and TeV scale gravitino requires low scale inflation [21]. 

This single-throat scenario is certainly the simplest possibility, and some variant of 
it might work for real world provided that supersymmetry generates further hierarchy 
from the inflation scale down to TeV scale. On the other hand, with a single throat 
scenario like this it is a little bit unclear how supersymmetry would be broken at the 
right scale and how to generate the small and positive cosmological constant which is 
observed in today's universe. For model building purpose, it gives us more room to 
consider flux compactification scenarios with more than one such throats. Here we 
comment on new issues in reheating such a multi-throat brane world. 

One immediate fact is that there is never a thermal equilibrium between any such 
pair of throats, due to the large redshift factors. Classically the localized degrees of 
freedom in any of the throat are effectively confined to the bottom of the throat, and 
can interact with those in another throat only via highly suppressed mixing operators. 
Wavefunction of such localized particles of mass m may have exponentially small 
tail toward top of the throat {m/MpY for some positive number d, so interaction 
amplitude between two different throats should carry suppression {m/MpY'^. In 
turn, the associated probability goes like {m/MpY^. Most crudely, the time scale of 
such interaction is then 



A hint on what the number d might be can be found in a work by Dimpopoulos 
et.al. 113 EO]; which estimated tunnelling rate from one Randall-Sundrum universe 
to another glued at the UV brane. The most optimistic estimate one finds there is 
d = l. 

We expect even lower tunnelling rate for angular modes. This extra suppression 
for angular modes can also be seen from the behavior of wavefunction outside the 
throats as well. In between throats, that is, in the "big" CY volume, we can put the 
warping to be almost constant, H{y) ~ Hq, and the wave equation reads as 



Thus in between throats solution of the KK wave equation (j2.19p is trivial, for L = 




(5.1) 



mKK\L — 



, 



(5.2) 
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modes it is simply a slow sinusoidal form (in fact more or less a small constant 
since itikkR is very small), while for angular KK modes it is exponential in 
^mKK;L ~ e^^^. As wavefunction traverses the classically forbidden region, it will be 
more strongly suppressed because of this behavior. For actual decay exponent with 
smooth Calabi-Yau (rather than singular boundary condition), we must do more 
careful analysis of the whole wavefunction. 

While the rate of decay of angular KK modes is seemingly thus more suppressed, 
this will probably lead to an overestimate of the lifetime of actual angular modes. 
One problem is that there is possibility that higher angular modes will decay a pair 
of lower ones classically without having to tunnel. To settle this, one must study 
lowest lying KK spectrum and dependence thereof. The other problem is that 
since the angular isometrics are approximate, the dominant decay channel might be 
to change the angular momentum to lower value within the same throat and then 
tunnel to another throat. For this, one must take into account perturbation of type 
that lead to ()4.18j) . Thus, one should not take this estimation of the net decay rate 
too seriously for larger values of L^. 

Let us return to the formula ()5.H) . With = 1, this time scale is clearly much 
larger that ~ Mp/m? at the exit of inflation, and thus the two throats are not 
in thermal equilibrium. While the 1/H increases with later evolution, so does At at 
much faster rate. Any exchange of energy between two such throats are possible only 
via possible decay of heavy particles in one throat into lighter ones in the other. 

This warns against one easy mistake in dealing with multi-throat cases. In early 
universe one often invokes equipartition principle and assume that energy is deposited 
in each and every degree of freedom in equal amount. A variant of this may be 
employed to argue that, if throat 1 comes with many more degrees of freedom than 
throat 2, then more energy is deposited to throat 1. However we are also familiar 
with the fact that some subsystem can be frozen away from thermal equilibrium and 
will evolve on its own if its interaction rate with the rest is much smaller that H . 
The above estimate of interaction between distinct throats tells us that such naive 
counting based on intuition from thermal equilibrium physics cannot be trusted. 

One can make this a little more precise using a simple quantum mechanics. Con- 
sider throat which drove inflation and contains massive KK modes of scale m after 
branes annihilation, and a longer throat 1 with mass scale -C m. The dominant 
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Decay 



Figure 3: KK modes in the inflation throat can decay to another throat via quantum 
oscillation. While the oscillation amplitude depends on mass eigenvalue distribution 
in the 2nd throat, it will be also suppressed further if the state in throat 1 has a large 
width. 

decay channel of heavy particle of mass m in throat is found by realizing that, since 
the two throats are not completely disconnected, KK modes of throat will mix and 
oscillate quantum mechanically with KK modes of throat 1 with similar mass, say 
m' ~ m ^ With a generic form of the mixing mass matrix, 

we expect e ~ m{'m/Mp)'^'^ <C m. The oscillation of the initial state in throat will 
induce amplitude of KK modes in throat 1 which is roughly 

(l|e-^^*|0) ~ e-^"^* ^ sin (Ami) (5.4) 

with Am ~ mi — m2- Part of the quantum state oscillating into the throat 1 will 
then decay since there are many lighter states in throat 1. Naively one may think 
that having many decay channel in throat 1 is very helpful in transferring energy into 
throat 1 from throat 0, since there are a lot of phase volume in throat 1. 
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However, to maintain consistency, the above mass matrix must be modified to 
reflect the decay width of the state in throat 1. 

j?=('" -r) (5.5) 

With number of states available as decay products in throat 1, the width should 
scale like F ~ Xm'N with a small but finite coupling constant A. With the new H 
that has the decay process built in, we can estimate how fast the initial state in throat 
lose its amplitude, 

e-'"'\0) ~ e-K-'"^' ( (1 -O(e'))|0) + ^sinAmill)^ (5.6) 

V r — iAm J 



where the effective width 



of the state |0) due to the oscillation is actually suppressed as the width F of state 
|1) increases. 

While somewhat counter-intuitive, this can be understood from the basics of the 
quantum oscillation. Two states mixes with each other well if they are almost de- 
generate while the mixing is suppressed by the mass difference. What this simple 
computation tells us is that the suppression depends on difference in the complex 
mass. Another way to see this is that the decay of state |0) occurs at second order 
in the perturbation theory. It must flrst oscillate into state |1) and experience the 
decay width F and then come back to |0). The standard 2nd order perturbation due 
to the mixing Hqi — e gives 

i^mmmi (6.8) 

Hii — Hoo 

as the width, which is exactly F. 

Finally we must take into account that there arc also roughly N number of states 
that can be used for such oscillation, so the net width of state |0) is more like 

A^f (5.9) 

Assuming large A^, so that Am is smaller than F, we find 

NT ~ —N ~ — (5.10) 
F Am 
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Figure 4: KK modes in the inflation throat deposit energy to lower energy throats. 
The branching ratio will be largely determined by how throats are distributed in the 
internal manifold with respect to the inflation throat and less sensitive to the field 
content of each throat. Energy deposited in throat 1 will later decay to throat 2, but 
at much more suppressed rate because its mass scale is far lower than that of the 
inflation throat. 

which shows that, to leading approximation, the effect of having a large number of 

state in the low energy throat completely washes out.-^- Any throat that is significantly 

longer that the inflation throat would receive energy as dictated by the mixing of 

typical KK wavef unctions, irrespective of how many decay channel each comes with. 

In particular, having a large number of open string decay channel in throat 1 will 

make matters worse, since these will contribute to F but not to the mixing. 

Such a universal nature of decay into different throats will make the reheating 

process quite delicate. For instance suppose the compact manifold involves three 

different throats. We will label them as the inflation throat, 1 the SUSY breaking 

throat, and 2 the standard model throat, with the mass scales m ^ //i ^ /X2, 

possibly with large number of additional open string degree of freedom in throat 2. 

Above consideration tells us that energy deposit into throat 1 and throat 2 are largely 

determined by each mixing mass matrix element e's and initial energy level m, and 

■^This estimate would not apply to relatively stable state of lower mass ~ in throat 1 at its 
bottom. However the corresponding e would be far more suppressed at ~ (m/x/Mj,)'', since the 
state |1) in that case would be living at the bottom of throat 1. 
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independent of /ii 2- This would result in significant deposit of energy in throat 1 as 
well as in throat 2, provided that the process is fast enough. In terms of ()5.1|) . d = 1 
with conventional inflation scale would suffice. 

With /ii at some intermediate scale, furthermore, whatever energy deposited there 
will behave as a massive particles as universe expands and cools down. Compared to 
the radiations in throat 2 (mostly photons and other light standard model degrees 
of freedom), the density of such massive particles dilute rather slowly. Even a small 
amount of deposit in throat 2 will quickly overtake and dominate energy density of 
universe rendering reasonable BBN impossible. While energy deposited in throat 1 
will eventually decay to throat 2, its decay width would be dictated by (fii/Mp)^'^ 
which is much smaller that {m/ MpY'^ . Relics in throat 1 would be very difficult to 
remove. We believe that this would cause a very serious problem whenever we build 
models based on KKLMMT-like inflation with warp factors playing crucial roles in 
hierarchy, determination of cosmological constant, and low scale of supersymmetry 
breaking. 

Notice that this problem with low-laying KK modes is independent from the 
angular KK modes. While problem with angular KK modes is there only if the base 
manifold of the throat happens to be quite symmetric, this problem is much more 
generic, provided that there are some intermediate-length throat containing a hidden 
sector. 

6 Summary 

It is fair to say that string theory inflation is at a crossroads. There are several 
potential directions for building realistic inflationary models. The end point of all 
inflationary models is reheating of the universe. Successful reheating means almost 
complete conversion of inflaton energy into thermalized radiation without any dan- 
gerous relics, in order to provide a thermal history of the universe compatible with 
Big Bang Nucleosynthesis (BBN), baryo/leptogenesis, and other observations. All of 
this provides us with constraints on string theory inflationary scenarios. 

In this paper we investigated reheating after brane-anti brane annihilation. The 
starting point was the KKLMMT model of warped brane inflation with a background 
KS geometry. Control of the calculations is possible in the regime of low curvature KS 
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geometry with all scales (like the radius at the bottom of the throat) larger than the 
string scale. Apart from the justification of computation, we do not have much reason 
to believe these constraints. In fact, for large to generate hierarchy, one may 

prefer i?_ nearer to the string scale fill. KS geometry also admits isometrics of the 
inner manifold. Another important parameter is the warping factor which defines 
the energy scale of inflation. The single throat scenario with low energy inflation just 
above the SM scale is the simplest possibility if there are alternative mechanisms to 
generate cosmological fluctuations (like modulated fluctuations). Otherwise we can 
have a multiple throats scenario with inflation and the SM sector located in different 
energy scale throats. Thus, the models bifurcate into single throat or multiple throat 
scenarios, with different ranges of parameters to stay in or out of low curvature 
geometry. 

We followed the way the energy of annihilated branes transfers into SM radiation. 
The first step after annihilation is the excitation of closed string loops, located near 
the bottom of the throat, where the local string mass scale is reduced by the warp 
factor e"^". They quickly decay into lOd gravitons (and particles in their supermul- 
tiplets). From a four dimensional perspective, lOd gravitons are manifested as usual 
4d gravitational radiation and massive KK modes. 

The generation of 4d gravitational radiation is a universal prediction of string 
theory brane inflation. In principle, it is a source of concern because BBN excludes 
more than few percent of radiation in forms other than photons and 3 light neutrino 
species. Massless radiation relics today are photons with the fraction of the total en- 
ergy density ^ 5 x 10^^ and very light neutrinos with Q^, ~ 1.6 x 10^^. Successful 
BBN potentially allows an extra few percent of energy density in light species, for in- 
stance background relic gravitons. Hence the energy density of relic gravitons cannot 
exceed the limit Qgw < 5 x 10^^. For instance, in brane inflation models without 
warping and just a pair of D3 — D3 branes, overproduction of gravitons is a problem, 
which can be cured by annihilation of D3 with a stack of D3 branes to dilute the 
abundance of gravitons by 1/A^. In warped geometry, the situation with gravitons is 
very different. The decay rate of closed loops into gravitons (i.e. homogeneous KK 
modes) is suppressed by Mp, while decay into KK modes is suppressed by the local 
string mass e'^^/^/a' at the bottom of the throat, which is much larger than Mp. 
As a result, the fraction of energy deposited into gravitons is reduced by the factor 
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e-22/0. For inflation at 10^^ GeV this is e-^^/o ^ iq-s. We therefore conclude that 
the warped brane inflation scenario predicts a residual background of relic gravita- 
tional radiation at the level ^gw ~ e~^^°. The wavelengths are located around the 
energy scale of closed loop excitations after inflation, redshifted by the expansion of 
the universe, and up to some numerical factors close to that of the relic photons, 
around mm scale. 

The bulk of the energy after decay of closed string loops is in KK modes. They 
interact and reach thermal equilibrium with a temperature comparable with the infla- 
tion energy scale. A much more serious concern comes from the angular KK modes. 
Indeed, the 6d compact manifold in the KS throat geometry has isometry, near the 
bottom X S'^ isometry. Conservation of angular momentum associated with these 
isometrics forbids complete decay of angular KK modes, but allows only their anni- 
hilation. As a result a residual amount of massive angular KK modes freezes out in 
an expanding universe, posing a problem for the scenario. 

In the context of CY compactification isometrics are only approximate, which al- 
lows the angular KK modes to decay. However, the width of their decay is suppressed 
by the factor [m/MpY . While this gives a fast enough decay time for angular KK 
modes in 10^^ GeV inflation, it is too long a time for angular KK modes in the SM 
throat. 

In fact, there are two independent KK mode problems. One has to do with 
the approximately conserved angular momentum which we just discussed. This is 
there for both the single and multiple throat cases as long as a long throat has an 
approximate isometry. The other problem comes from a low-lying KK mode created in 
some intermediate throat (low-scale inflation throat or susy-breaking hidden-sector 
throat) in multi-throat scenarios. Their decay is via tunnelling only, and is thus 
suppressed. This suppression has little to do with any isometry and is thus more 
problematic, at least for generic multi-throat cases. 

It is an interesting exercise to calculate the cosmological abundance of residual 
KK modes, which we considered only briefly for the single throat case in Section 4. 
Let us mention the cosmological relevance of the issue. The energy density of massive 
relics ^kk dilutes slower than that of radiation, therefore ^kk sooner or later begins 
to dominate the expansion of the universe. Special conditions are needed to tune 
^KK to be the Cold Dark Matter component, to constitute today Qcdm — 0.27. If 
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KK particles are abundant in the very early universe it becomes matter dominated 
too early to allow successful BBN. This is a typical situation with SUSY moduli 
fields and gravitinos. A similar problem is inherited by simple versions of rolling 
tachyon inflation |3] . On the other hand, it is fair to keep a possibility that Qkk is a 
candidate for the CDM, if KK particles are long-lived and are created in exactly the 
right amount ;51] . 
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